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■ In recent years an increasing attention is paid to the so-called correlated disorder in low- dimensional 



disordered systems. The interest to this subject is mainly due to two reasons. First, it was found 
O : that specific correlations in a disordered potential can result in a quite unexpected anomalous prop- 
, erties of scattering. Second, it was shown that such correlations can be relatively easily constructed 
^ ' experimentally, at least, in the one-dimensional Anderson model and in Kronig- Penney models of var- 
i/-^ ious types. Therefore, it seems to be feasible to fabricate random structures with desired scattering 
O ■ properties, in particular, when one needs to suppress or enhance the localization in given frequency 
windows for scattering electrons or electromagnetic waves. In addition, it was understood that in 
many real systems the correlated disorder is an intrinsic property of underlying structures. One of 
the most important examples is a DNA chain for which strong correlations in the potential have 
shown to manifest themselves in an anomalous conductance. Thus, the subject of correlated disorder 
is important both from the theoretical viewpoint, and for various applications in physics. 

\ The key point of the theory of correlated disorder is that the localization length for eigenstates 

^ ■ in one-dimensional models absorbs the main effect of correlations in disordered potentials. This fact 
^ ■ was known since the earliest analytical studies of transport in continuous random potentials, however, 
until recently the main interest was in delta-correlated potentials, or in potentials with a Gaussian- 
type of correlations. On the other hand, it was shown [U [21 El S] that the most interesting effect 
is related to specific long-range correlations that can be fabricated in practice. In particular, it was 
demonstrated that in the standard one-dimensional Anderson model one can observe effective mobility 
edges in the energy spectrum, when the pair correlator computed along the disorder is of specific form. 
The important feature of the mobility edge Uc in the frequency spectra of traveling waves, is that 
it separates the region of strongly localized states from that of very extended states. This property 
is important for the construction of potentials with selective transmission or refiection. From the 
experimental viewpoint, many of the obtained results may have a strong impact for the creation of a 
new class of electron nanodevices, optic fibers, acoustic and electromagnetic waveguides with selective 
transport properties. 
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In spite of an asymptotic character of theoretical results obtained for infinitely large samples 
and weak disorder, the analytical predictions were found to work relatively good for a strong disorder 
as well. The first experimental study [21 El C] of both the suppression and enhancement of localization 
due to correlated disorder, was performed on single-mode electromagnetic waveguides. It was shown 
that in the case of statistically correlated point-like surface scatterers one can create the controlled 
frequency windows of the enhanced transmission, or windows with a very strong reflection. The 
important point is that in spite of many experimental imperfections, such as very strong absorbtion 
or small number of scatterers with large amplitudes, the selective transport was clearly observed in 
accordance with the theory. 

For single- mode waveguides the problem of the surface scattering can be reduced to that for one- 
dimensional disordered models. For this reason the methods and results obtained for the latter case can 
be directly applied for the waveguides [H |9l [TOl [TTl [12] . The situation is principally different for many- 
mode waveguides. The problem of wave propagation through such systems with corrugated surfaces 
has a long history and till now remains a hot topic in the literature. This problem naturally arises in 
the analysis of spectral and transport properties of optics fibers, acoustic and radio waveguides, remote 
sensing, shallow water waves, multilayered systems and photonic lattices, etc. [I3l [TH UHl HEl [HI HE] . 
Similar problems emerge in quantum physics, when describing the propagation of quasi-particles in 
thin metal films and semiconductor nanostructures, such as nanowires and strips, superlattices and 
quantum-well-systems [l9l[20l[2ll[22l[23l[2ll[25l[26l[27]. Recently, new theoretical results [281 129] 
have been obtained for one-dimensional disordered models describing photonic crystals, bi-layered 
metamaterials and electronic superlattices. The experiments performed on microwave guiding systems 
with intentionally randomized model parameters [30], have confirmed the predictions of the theory. 

As is well established, the scattering from corrugated surfaces results in the diffusive transport 
[3T1 l32l [331 EH [351 EHl [371 ES [39], as well as in the effects of strong electron/wave localization 
[ini [HI [121 [131 [m [ISl [ini [13 \M\ ■ correspondingly, the eigenstates of periodic systems with corrugated 
surfaces turn out to have a chaotic structure [19]. Recent numerical studies of quasi-one- dimensional 
surface- disordered systems [IHl [13 [H] have revealed a principal difference from those known in the 
standard models with bulk random potentials [23] . Specifically, it was found that transport properties 
of quasi-one-dimensional waveguides with rough surfaces essentially depend on many characteristic 
lengths. In comparison, for the bulk scattering all transport characteristics depend on one parameter 
only, which is the ratio of the localization length to the size of a sample (the so-called single-paramter 
scaling) . 

The situation in many-mode waveguides in the presence of long-range correlations was found 
to be quite tricky [lOl [501 EI]- It was shown that the long-range correlations, on the one hand, give 
rise to a suppression of the interaction between different propagating waveguide modes. On the other 
hand, the same correlations can provide a perfect transparency of each independent channel, similar 
to what happens in the one-dimensional geometry. The number of independent transparent modes 
is governed by the correlation length and can be equal to the total number of propagating modes. 
Therefore, the transmission through such waveguides can be significantly enhanced in comparison 
with the case of uncorrelated surface roughness. 

It should be stressed that the main results in the theory of surface scattering were obtained 
for random surfaces with fast-decaying correlations along the structures. Therefore, it is of great 
importance to explore the role of specific long-range correlations in surface profiles, using the results 
found for one-dimensional systems with correlated disorder. Apart from the theoretical interest, 
this problem can be studied experimentally since the existing experimental technics allow for the 
construction of systems with sophisticated surfaces [52l [53] . 

We would like to note that in order to focus on the role of long-range correlations in surface 
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profiles, in this review we do not discuss the so-called square- gradient mechanism of scattering. As 
was recently shown, this mechanism emerges due to a quite specific dependence of the back scattering 
length on the second derivative of scattering profiles. The theoretical aspects of the square-gradient 
scattering and its possible experimental implications can be found in Refs. |54] . 



2 Surface- corrugated waveguide 

As a physically plausible and commonly used model to study the multiple surface scattering, we 
consider an open plane waveguide of length L and average width d with perfectly conducting lateral 
walls. It is naturally to imply the waveguide length be much greater than its width, d <^ L. Such a 
system is called quasi-one- dimensional one. The x-axis is stretched along the structure and the z-axis 
is directed in transverse direction. The lower boundary of the waveguide is supposed to have rough 
(corrugated) profile z = ^(x), slightly deviated from its flat average z = 0. The upper boundary is 
taken to be fiat, z = d, (see Figured]). Thus, the surface-corrugated guiding system occupies the area 
defined by the relations, 

-L/2<x< L/2, ^(x) ^z^d. (1) 
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Figure 1: (Color online) Plane waveguide with lower corrugated edge. 

The random function ^ (x) describing the surface roughness is assumed to be statistically homo- 
geneous with the following characteristics, 

(e(x)) = 0, (eix)) = a^ iam^')) = cr'Mx - x'). (2) 

The angular brackets (. . .) stand for ensemble averaging over the disorder, i.e. over different real- 
izations of the random surface profile ^{x), or for the spatial average over the coordinate x of any 
prescribed realization. Both types of the average are assumed to be equivalent due to ergodicity. The 
variance of ^(x) is denoted by and, consequently, a is the root-mean-square roughness height. The 
binary (two-point) correlator W(x) is normalized to its maximal value, W(0) = 1, and is assumed to 
decrease with increasing |x| on characteristic scale termed the correlation length. 

In what follows we consider weak surface scattering for which the corrugations are small, a <^ 
d. This limitation is common in the surface scattering theories based on appropriate perturbative 
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approaches [13]. As is known, for weak scattering all transport properties are entirely determined by 

the roughness power spectrum W{kx), 



-^exp{tk,x) W{K), (3a) 

-oo ^TT 

/oo 
dxexp{-ihx)W{x). (3b) 
-oo 

Since the pair correlator W{x) is a real and even function of the coordinate x, its Fourier transform 
W{kx) is even, real and non-negative function of the lengthwise wave number k^. Note that the 
condition W(0) = 1 is equivalent to the following normalization for W{kx), 



J^W{kx) = l. (4) 

Since in the x-direction the system is open, at x = ±L/2 we assume the radiative boundary 
conditions. In the transverse 2;-direction the zero Dirichlet boundary conditions are applied to both 
lateral walls, z = ^{x) and z = d. Thus, the analysis of the surface scattering in our model is reduced 
to the study of the following two-dimensional boundary-value problem. 



^(a;,z = ^(x)) = 0, m{x,z = d) = Q. (5b) 

Here k = u/cis the total wave number for electromagnetic wave of frequency u and TE polarization, 
propagating through a waveguide with perfectly conducting walls. Note that in contrast with the 
bulk-scattering, here the wave equation does not contain any scattering potential since the source of 
the scattering is the roughness of boundaries. 

In ideal waveguide with fiat walls, ^(x) = 0, the solution of the problem ([5]) has the canonical 
form of normal waveguide modes. 



1 /"7inz\ 

W„,±(x, z) = —= sm — — exp(±zfc„x). (6) 
\ -Kd \ d ) 



Txd 

Here the integer n = 1, 2, 3, . . . enumerates the normal modes ([6]) with the transverse wave number 
^2 = Tin/d. The lengthwise wave number k^ = ±kn for the nth mode is. 



fc„ = y/k^ - {nn/dy. (7) 

Evidently, the transport properties depend only on normal modes that can propagate along the 
waveguide, i.e. have real value of kn- As follows from Eq. ([7]), such propagating modes have indices 
n ^ Nd, and their total number is equal to the integer part [...] of the ratio kd/n, 



Nd = [kd/Ti]. (8) 

The waveguide modes with indices n > have purely imaginary wave numbers kn- These evanescent 
modes decay exponentially fast on the scale of the order of wavelength. As one can see, the unper- 
turbed (fiat) waveguide is equivalent to the set of A^^ one- dimensional non-interacting conducting 
channels occupied by the corresponding propagating modes. 
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3 Single-mode structure 



Keeping in mind the relevance of wave scattering to the Anderson locahzation, let us, first, consider 
a single-mode waveguide. In this case the mode parameter kd/ir is restricted by the relation 1 < 
kd/TT < 2, and the number of conducting channels equals one, A^^^ = 1. The transmission through 
such waveguide depends on the lengthwise wave number ki, 



k, = VP - {ix/df. (9) 

All other waveguide modes with n ^ 2 are evanescent and do not contribute to the transport prop- 
erties. From the single-mode condition it follows that the wave number ki is confined within the 
interval, 

< kidl-K < V?>. (10) 
Note that the weak surface-scattering condition a d, leads to the inequality kicr ^ 1. 

As was shown in Refs. |12], the transport problem (|5]) for the surface-disordered single- mode 
waveguide is equivalent to a one- dimensional disordered model, 

' d^ 



— + k'- V{x 



ip{x) = (11) 



where ki stands in place of k. In such a description, the effective potential V{x) has the form, 

V^(^) = I ax), (12) 

which is entirely determined by the rough surface profile ^(x). 

As one can see, the surface scattering in one-mode waveguides is equivalent to the bulk scattering 
emerging in one-dimensional disordered systems. The latter problem can be solved with the use of 
well-developed methods, such as the perturbative diagrammatic technique of Berezinski [20j, the 
invariant imbedding method |42l |55l [56] or the two-scale approach [44l |57] . All these methods allow 
one to take adequately into account the effects of coherent multiple scattering from the corrugated 
surface giving rise to the Anderson localization. 

The main theoretical result is that the average transmittance (T) as well as all its moments 
(T*), are described by the universal function. 



{T'{L/Lioc)) = exp ( ^] 

2a/ vr \2LiocJ \ ILiocj 

s = 0,±1,±2,... . (13) 

This function depends solely on the ratio L/Lioc between the waveguide length L and the localization 
length Lioc (see, for example, Ref. j21]). 

In one-dimensional disordered systems a wave can be scattered either forward or backward. 
However, it was shown that the transport properties are determined exclusively by the backscattering 



5 



while the forward scattering has no influence. Therefore, the quantity Lioc is, in fact, the hacks cattering 
length emerging in an infinite one-dimensional structure. It is important that the inverse value L;~J 
can be associated with the Lyapunov exponent appearing in various transfer matrix approaches |2T] . 
In the latter description the Lyapunov exponent gives the average rate of decrease of the wave function 
{ip{x)) away from the center of its localization. 

It should be stressed that the dependence of the transport properties on the ratio L/Lioc mani- 
fests a principal concept of one-parameter scaling that constitutes the phenomenon of one-dimensional 
Anderson localization. The non-trivial point of this concept is that in order to describe the transport 
properties of finite samples of size L, it is sufficient to know how the wave function is localized in an 
infinite sample with the same disorder. 

From Eq. f|T3|) one can derive relatively easy the expressions for low moments of the transmit- 
tance T. Specifically, at s = 1 one obtains the average transmittance {T{L/Lioc)) , 



{T{L/Lioc)) = exp ( 

2V7r \2LiocJ V '^LiocJ 

The second moment, s = 2, is important to get the variance of the transmittance. It can be shown 
that for Ljoc ^ ^i the variance is of the order of squared average transmittance itself. This means 
that for strong localization the transmittance is not self-averaged quantity. Hence, by changing the 
length L of the waveguide, or the disorder itself, one should expect very large fiuctuations of the 
transmittance. Such fiuctuations are known as the mesoscopic fluctuations that are characteristic of 
strong interference effects on macroscopic scale. 

In order to properly characterize the transport properties of one- dimensional structures for any 
degree of localization (ratio L/Lioc), one should refer to the self- average logarithm of the transmittance, 

(lnT(L/L;„,)) = -2L/Lioc- (15) 

This result is consistent with an exponential decrease of the transmittance averaged over the so-called 
representative (non-resonant) realizations of the random disorder |21] . 

{T{L/Lioc))rep = exp{-2L/Lioc)- (16) 

Note that Eq. f|T5l) quite often is used as the definition of the localization length Lioc itself. It 
is highly non-trivial that by exploring the transmission properties of finite samples, one can extract 
the localization length that is defined for infinite samples. This fact is again the manifestation of the 
one-parameter scaling. 

In accordance with the scaling concept, in the one- dimensional geometry there are only two 
characteristic regimes corresponding to the ballistic and localized transport. 

(i) The ballistic transport occurs if the localization length Lioc is much larger than the system 
length L. In this case one- dimensional structure is practically fully transparent since its average 
transmittance is close to one, 

{T {L/ Lioc)) ^l-2L/ Lioc for L,,, > L. (17) 
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This asymptotic expression results from both expressions, (IT^ and (flGll . 

(ii) Otherwise, the disordered structures exhibit the localized transport, when the locahzation 
length Lioc is smaller than the sample length L. In this case the average transmittance (HM is 
exponentially small, 



for Lioc < L. (18) 

This means that in the localization regime the disordered one-mode waveguides perfectly (with an 
exponential accuracy) reflect the incoming waves. 

According to the above expression the transmission exponentially decreases on the scale L ^ 
2Lioci with an additional power prefactor. In contrast, the transmittance ( IT6ll has an exponential 
dependence with much faster decrease on the scale L ^ Lioc/ 2. This fact can be explained as follows. 
The main contribution to the asymptotic form f lTSj) for the average transmittance ( !T^ is given by 
resonant realizations of the random potential V{x). For these realizations the transmittance is almost 
equal to one, however, they have an exponentially small probability. On the other hand, for represen- 
tative realizations (most probable, but non- resonant) the transmittance is described by formula ([TB|1 . 
This effect is peculiar to the mesoscopic nature of Anderson localization. 

The expressions ( lT3ll - (JT8l) are universal and applicable for any one-dimensional system with 
a weak static disorder. As one can see, in order to describe transport properties of finite structure, 
one needs to know the localization length Lioc- According to different approaches |20l |2T1 |55l |56l |57] . 
the inverse localization length for any kind of weak disorder is determined by the 2/c-harmonic in the 
randomness power spectrum S{kx) of the scattering potential V{x), 



L~lik) = Si2k)/8k'; (19a) 
{V{x)V{x')) =C{x-x'), (19b) 

/oo 
dx exp{—ikxx) C (x) . (19c) 
-oo 

For elastic backward scattering the wave vector k conserves its value k, changing the sign, | AA;| = 2k. 
Accordingly, Eq. f lT9|) reflects the fact that the localization length is defined by the backscattering 
length only. 

The expression ( !T9|) indicates that the global properties of the wave transmission through one- 
dimensional disordered media depend on the two-point correlations in the random scattering potential. 
Therefore, if the power spectrum S{2k) is very small or vanishes within some interval of the wave 
number k, then the localization length Lioc appears to be very large {Lioc ^ L) or even diverges. 
Evidently, the localization effects can be neglected in this case, and the structure, even of a large 
length, is fully transparent. This means that, in principle, by a proper choice of the disorder one can 
design the disordered structures with selective [anomalous) ballistic transport within a prescribed 
range of k. 

Taking into account the form (!T2|) of the potential and its correlation properties ([2]), from 
Eqs. ( !T9|) one can readily derive the following explicit formula for the localization length in the single- 
mode waveguide [42j , 
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Since the potential fll2l) is entirely determined by the rough surface profile ^(x), the localization 
length fl20|) is specified by the roughness power spectrum W{kx)- Therefore, by a proper fabrication 
of a random profile ^ (x) with specific long-range correlations, one can arrange a desirable anomalous 
transport within a given window of k = u/c inside the single- mode region ffTOj) . 

Before we start with a practical implementation of the expression (!20|) for the inverse localization 
length, it is worthwhile to clear up some points. First, one should stress that this result, as well as 
Eq. (fT9l) obtained by different methods, is an asymptotic one. This means that the higher terms are 
non-controlled; however, they can be neglected in the limit cr^ — )■ 0. Second, the main assumptions 
used in the derivation of Eq. (I2UI) are based on the validity of averaging over different realizations 
of disorder (see, for instance, Ref. [57]). The condition for such an average is that two lengths, Lioc 
and L, are much larger than two other characteristic lengths, the wavelength fcf ^ and the correlation 
length determining the maximal value of the power spectrum W{2ki). One should stress that for 
any finite values of ki and kc, this condition can be always fulfilled due to an asymptotic character of 
Eq. fl2U]) (or, the same, due to small value of a^). 



4 Design of random surface profile with predefined correla- 
tions: Convolution method 

From the above consideration it is seen that, in principle, by a proper choice of surface disorder one 
can artificially create the systems with selective transparency or reflectivity. Thus, the important 
practical problem arises of how to construct a corrugated surface profile from a predefined roughness 
power spectrum W{kx). This problem can be solved by employing a widely used convolution method 
that was originally proposed in Ref. [58]. The modern applications of this method for generation of 
random structures with specific correlations, including long-range non-exponential correlations, can 
be found in Refs. [591 ESI EOl EU |62l [63], [61] and in other papers cited in this chapter. 

The method consists of the following steps. First, having a desirable form for the power spectrum 
W{kx), we derive the modulation function (3(x) whose Fourier transform is W^^'^^k^), 

/°° dk 
-^exp(^fc,.x) W^'\kx). (21) 
-oo 27r 

Then the random surface profile ^{x) is generated as a convolution of a white-noise Z{x) with the 
modulation function /3(x), 

/oo 
dx' Z{x-x')(5{x'). (22) 
-oo 

The delta-correlated random process Z{x) is determined by the standard relations, 

(Z(x)) = 0, {Z{x)Z{x')) = 5{x~x'), (23) 

and can be numerically created with the use of random-number generators. Here, 5{x) is the Dirac 
delta-function. 
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The expressions (122|) and (!23|) give the solution of the inverse scattering problem of constructing 
random roughness from its power spectrum. Note that this method is valid in the case of a weak 
disorder only. That is why only the binary correlator is involved in the construction of ^ (x) while the 
higher correlators do not contribute. Note also that the profile obtained by the proposed method is 
not unique. Indeed, there is an infinite number of realizations of delta-correlated noise Z{x) that give 
rise to different profiles ^{x) having the same power spectrum W{kx). 

The importance of this method is due to a possibility to obtain profiles resulting in a sharp 
transition between ballistic and localized transport, when changing the wave number k. In this case 
the corresponding power spectrum W{kx) abruptly vanishes at prescribed values of kc- This means 
that the binary correlator >V(x) has to be a slowly decaying function of the distance In other 
words, the corresponding corrugated surface profiles ^(x) should be of specific form revealing long- 
range correlations along the structure. Because of an abrupt character of transmission properties, the 
transition point kc can be regarded as an effective transparency edge. 

As is pointed out above, a statistical treatment is meaningful if the scale of decrease k~^ of 
the correlator W{x) is much less than both the sample length L and localization length Lioc- In this 
connection one should stress that long-range correlations we speak about, do not assume large values 
for the correlation length k~^. Indeed, the simplest correlator, >V(x) = sin kcx/kcX has finite scale 
k~^ of its decrease and k~^ can be quite small, see examples below. On the other hand, an effective 
width of the transparency edge is determined by the product, {Liockc)~^, not by k~^, and turns out 
to be very small (for details, see, e.g., Ref. ^3). One can say that the sharpness of the transition is 
defined by the form of the pair correlator rather than by the value of its correlation length. 

Note that the systems with very complicated scattering potentials are not exotic. For example, 
bulk random potentials have been constructed in the experiments [3 O [6l [53] , while rough surfaces 
with rectangular power spectrum have been fabricated in the study of the backscattering enhancement 



5 Gaussian correlations 

As was mentioned above, the general expression fl2U]) for the localization length Liodki) indicates that 
all features of the wave transmission through the surface-disordered single-mode waveguide depend 
on two-point correlations in the surface profile. In order to demonstrate how to realize the properties 
of long-range correlated disorder, let us first consider the surface roughness ^(x) with a widely used 
Gaussian correlator. 



This correlator exponentially decreases on the scale of the correlation length k^ ^. 

Using the convolution method ( 12T]) - ( l23l) one can obtain that the corrugated surface profile 
^(x) with the correlation properties ( 124|) is described by the function. 



W{k,) 



>V(x) 



exp (— fc^x^) , 
V^k;'exp{-ej4kl). 



(24a) 
(24b) 




(25) 



Correspondingly, the inverse localization length ( 120|) takes the following explicit form. 
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^^-^^'>-^^Ad) {2k,Y ■ ^^^^ 



One can see that within the single-mode interval (fTOj) the localization length increases exponentially 
with ki from zero at fci = to a large value at ki = vrv^/rf- Clearly, in the vicinity of fci = 
the localization length is Lioc{ki) is much less than L and the waveguide is non-transparent. Thus, 
the localization regime ( |T8|) occurs within the whole single-mode interval, provided by the condition 

LiociT^V^/d) L at ki = 7r\/3/d. 

On the contrary, when Lioc{TT\/3/d) ^ L, one can observe the crossover from the localized 
transport fllSp to the ballistic one fll7l) . For the Gaussian correlations ([MD, both the crossing point 
where Lioc{ki) = L, and the crossover width depend on the values of k~^ and L. Because of a smooth 
character of the crossover, the crossing point cannot be regarded as the transparency edge. However, 
one can see that the longer the correlation length k~^, the larger localization length Lioc{ki). Hence, 
the larger is the ballistic region and more narrow is the crossover. One can conclude that, in general, 
the correlations suppress the localization. 

The surface profile ^(x) with Gaussian correlations admits the uncorrelated roughness of the 
white-noise type. Indeed, since cr'^k'^ = const, from Eqs. for k~^ — )■ one obtains the delta-like 
correlator and constant power spectrum. 



W^n{x) = V^k;^6{x), (27a) 
W^n{k.) = ^k-\ (27b) 

The convolution method (l^T]) - fl23p results in the following expression for the surface profile 

ax): 

e^„(x) = ^^Z(x). (28) 



According to Eq. ( 120|) . the localization length reads. 



The expressions ( |271) - (1291) can be considered as the asymptotical ones of the respective equations 
(El - ([26]) when (fci/fc^)^ < 1- 

The comparison of Eqs. f l29|) and fl26|) leads to the conclusion: the best way to observe localized 
transport is to employ an uncorrelated disordered surface. Indeed, the condition Lioc{ki) ^ L is 
stronger for the Gaussian correlations than for the case of white-noise profiles. On the other hand, 
for Gaussian correlations with a small value of kc, the ballistic regime ( !T7|1 can be realized even for 
such lengths L and wave numbers ki for which strong localization, ( !T8|) . takes place for delta-like 
correlations. Again, this fact confirms that the correlations suppress the localization. 



6 Two complementary examples of selective transparency 

As we discussed in Section IH the rough surfaces with prescribed two-point correlations can be con- 
structed with the use of the convolution method. Below we demonstrate the construction of surface- 
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disordered structures with selective transparency, by considering two examples of long-range correla- 
tions. 

(a) Let us first consider the waveguide which is non-transparent when the wave number ki is 
less than some value kc, and completely transparent for ki > k^. Such a behavior can be observed 
if the transition point (transparency edge) ki = k^ is located inside the allowed single-mode interval 
(COD, 



< kj/n < Vs. (30) 
For this case one can get the following expressions for the binary correlator >V(x) and power spectrum 

w{k^y. 

Wa{k,) = ^Q{2k,-\k,\). (31b) 



2k 



c 



Here Q{x) is the Heaviside unit-step function, Q{x < 0) = and 0(x > 0) = 1, and k^ is the 
correlation parameter to be specified. 

According to the recipe (12T]) - (!23ll . the surface profile with the above properties is given by the 
expression 

c(\ ^ r ^ ,ry, ,, sin(2fcexO 

Ux) = —== dxZ{x-x) . (32) 

Correspondingly, the inverse localization length has the step- down form, 

^'»<*')-i(i)'^- 

In line with this expression, as ki increases, the localization length Lioc{ki) also smoothly increases 
and then diverges at ki = kc- Thus, within the region < ki < kc the average transmittance 
{T{L/Lioc)) is expected to be exponentially small ( JTSll due to strong localization. The condition for 
strong localization to the left from ki = kc reads, 

>1. (34) 



Lioc{kc-0) Ankl \d 

Otherwise, inside the interval kc < ki < -n^/S/d a ballistic regime occurs with perfect transparency, 
{T{L/Lioc)) = 1. 

(b) The second example refers to a complementary situation when for ki < kc the waveguide 
is perfectly transparent and for ki > kc is non-transparent. The corresponding expressions for W(x) 
and W{kx) are given by, 

W,{x) = vr5(2M - ^^^^7^, (35a) 
W,{k,) = ^e{\k,\-2kc). (35b) 
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In this case the corrugated surface is described by a superposition of a white noise and the roughness 
of the first type, 



x' 



(36) 



Correspondingly, the inverse locahzation length is expressed by the step-up function. 



As a consequence, in contrast with the first case, here the surface-scattering localization length Lioc{ki) 
diverges below the transparency edge ki = k^. At this point Lioc{ki) sharply falls down to a finite 
value Lioc{kc + 0), and then smoothly increases with further increase of ki. In order to observe the 
localization within the whole region kc < ki < ix^^jd, one should assume that strong localization is 
retained at upper point fci = ix^jd of the single- mode region ffTOj) . 

L a^L /7r\4 ^ ^ 

> 1. (38) 



LiociirV^/d) l2TTk^ \d. 

Therefore, in this example the ballistic transport is abruptly replaced by strong localization at the 
transparency edge, ki = kc- 

One should stress that the surface profiles fl32p and fl36p with respective binary correlators and 
power spectra flHT]) and fl5^ are substantially different from the delta-correlated white noise f l25]) . fl27|) . 
and from random processes (l25|l with exponentially decaying Gaussian correlations i^^. Specifically, 
here the profiles are characterized by the variation scale {2kc)~^ and have the long power-decaying 
tales in the expressions for their two-point correlator. Such tails are originated from the stepwise 
discontinuity at the points k^ = ±2kc. The location of the transparency edge is defined by these 
discontinuity points and do not depend on other parameters, in contrast to the case of Gaussian or 
delta-like correlations. 

Now we demonstrate the above predictions by a direct numerical simulation. For this, the inverse 
localization length L^"^ was computed with the use of the Hamiltonian map approach developed in 
Refs. [H [2]. First, the continuous scattering potential ( [T2|) was approximated by the sum of delta 
kicks with the spacing 6 chosen much smaller than any physical length scale in the model. After, a 
discrete analogue of the ID wave equation ( ITTl) was analyzed numerically, with the lengthwise wave 
number ki in place of k and with the surface scattering potential V{x) given by Eq. (JT2l) . In this way 
the wave equation was expressed in the form of a two-dimensional Hamiltonian map describing the 
dynamics of a classical linear oscillator under the parametric noise determined by V{x). As a result, 
the analysis of the localization length was reduced to the computation of the Lyapunov exponent L^~J 
associated with this map (see details in Refs. [I1I2])- 

Numerical data reported in Figure |2l represent the dependence of the dimensionless Lyapunov 
exponent A = cqLJ'^]. on the normalized wave number K = ki/kcin the range < K < 2 corresponding 
to the single-mode interval. The normalization coefficient cq was set to have A = for the 

delta-correlated potential. Two surface profiles ^(x) were generated according to discrete versions 
of Eqs. f p2|) and fl36|) determining complementary stepwise dependencies of the localization length 
Lioc{ki) in accordance with Eqs. ( 133|) and f l37|) . 

One can clearly see a non-trivial dependence of A on the wave vector A', which is due to specific 
long-range correlations in ^(x). The data display sharp dependencies of A on A' when crossing the 
point A' = 1. Thus, by taking the size L of the scattering region according to the requirements ([3i 
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0.0 0.5 1.0 1.5 2-0 0-0 0.5 1.0 1.5 2.0 

K K 

Figure 2: Selective dependence of the rescaled Lyapunov exponent on wave number for two realizations 
of random surface with specific long-range correlations. To show the main effect of correlations, the 
complementary dependence of the Lyapunov exponent on K is shown: (a) eigenstates are localized 
for K < 1 and delocalized for K > 1; (b) complete delocalization for K < 1 alternates by strong 
localization for K > 1 (after [1]). 



or (1551) . one can arrange anomalous transport in the single- mode guiding structure as is predicted by 
the analytical theory. 



7 Random narrow-band reflector 

Let us now consider one more example of the binary correlator giving rise to the power spectrum of 
the following rectangular form, 



W{K) = ^ ^ [e(2A:+ - Ifc.l) - 9(2fc_ - |fc.|)] , (39b) 

< A;_ < A;+ < nVs/d. 

In these relations the factor 1/2(A;+ — A;_) stands to provide the normalization requirement ([1]), or 
the same, W(0) = 1. Such a power spectrum has been employed to create specific rough surfaces 
in the experimental study of an enhanced backscattering [52] . This spectrum was also used in the 
theoretical analysis of light scattering from amplifying media, as well as in the study of localization 
of plasmon polaritons on random surfaces 

In accordance with the convolution method fl2T|) - fl23|) . the random surface profile ^(x) with 
the above correlation properties can be obtain from the following expression: 

^/ \ f°° , f sin(2A;+x') - sin(2A;_a:') 



The peculiarity of such surfaces is that they have two characteristic scales, {2k+)~^ and (2A;_)~^. 
Consequently, the binary correlator and its power spectrum fl39|) are specified by two correlation 
parameters. 

From Eqs. (!20|) and (!39|) one can find the inverse localization length. 
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_ a' /vrx 6 e(fc^ - fc09(fci - fc_) 
^'-^^^^-7r(/t^.-A:_) {2k^^ " ^^^^ 

As one can see, there are two transparency edges at the points ki = A;_ and ki = kj^. The locahzation 
length Liociki) diverges below the first point, ki = k-, and above the second one, ki = k+. Between 
these points, for fc_ < fci < the localization length (HT]) has finite value and smoothly increases 
with an increase of wave number ki. 

Let us now choose the parameters for which the regime of a strong localization occurs at the 
upper transition point ki = where Lioc{ki) gets its maximal value. This automatically provides 
strong localization within the whole interval k_ < ki < k^. The condition to have this situation reads 

>1. (42) 



L^„,(fc+-0) A'Kkl{k+-k_)\d, 

As a result, there are two regions of perfect transparency for the waveguides of finite length L with the 
specified above surface profile. Between these regions the average transmittance (T) is exponentially 
small according to the expression (ITSl) . Due to this fact, the system exhibits the localized transport 
within the interval k_ < ki < k^ and the ballistic regime with (T) = 1 outside this interval. In the 
experiment one can observe that with an increase of the wave number ki, the perfect transparency 
below ki = k_ abruptly alternates with a complete reflection, and recovers at ki = From 
equations fHTj) and fH2l) one can see that the smaller value /c_|_ — /c_ of the reflecting region the smaller 
the surface-scattering localization length Lioc{ki) and, consequently, the stronger is the localization 
within this region. This remarkable fact may find important applications in creating a new class of 
random narrow-band filters or reflectors. 



8 Multi-mode waveguide 

Now we examine the correlated surface scattering in multi-mode waveguides, i.e. in waveguides with 
a large number > 1 of conducting channels, see expression ([8]). According to Landauer's concept 
[2l] . the total average transmittance (T) of any disordered quasi-one-dimensional guiding structure 
can be expressed as a sum of partial average transmittances {Tn) that describe the transport for every 
n-th propagating normal mode, 

{T) = J2(Tn). (43) 

n=l 

When all conducting channels are open, i.e. all T„ = 1, the total transmittance 0431) attains its 
maximal value equal to the total number A^^ of propagating modes. Therefore, our definition of the 
transmittance differs from the canonical one in which the maximal value of total transmittance is equal 
to one. Nevertheless, we shall use the definition (B3l) in order to clearly discriminate the intervals of 
the mode parameter /cd/vr with different numbers of the conducting channels and more clearly display 
the role of correlated surface disorder. In order to pass to the definition used in wave theories, one 
should simply divide our expression (H3ll over A^^. 

From general theory of quasi-one-dimensional scattering systems it follows that the transmission 
properties of any nth conducting channel (1 ^ n ^ A^^) are determined by two attenuation lengths, the 
forward scattering length Ln^ and hacks cattering length Ln^. For multi-mode quasi-one-dimensional 
waveguides with surface disorder the inverse scattering lengths are given by. 
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1 _,(^-(^,^^ 

-i^n " n'=l 

Here the lengthwise wave number /c„ is defined by Eq. (j7]). The results (jHj) and (H5|) can be obtained 
for the boundary- value problem ([5]) by the diagrammatic Green's function method [13], as well as by 
the technique developed in Ref. [40j. Also, these expressions can be derived by using the invariant 
imbedding method extended to quasi-one- dimensional structures [IE]. Note that in a single-mode 
waveguide with Nd = 1 the sum over n' contains only one term with n' = n = 1. Therefore, in 
this case the backscattering length hf^ is exactly equal to the single-mode localization length fl20|) . 



Lf^ = Lioc{h). see Eq. ([H. 

The expressions (jH]) and fH5|) manifest that, in general, both attenuation lengths are con- 
tributed by scattering of a given n-th propagating mode into all other modes. This is the case when, 
for example, a rough surface profile is either delta-correlated random process with constant power 
spectrum, W{kx) = const, or has fast decreasing binary correlator yV{x) and, correspondingly, slowly 
decreasing roughness power spectrum W{kx)- 

Another peculiarity is that the expressions (l44l) and (l45ll display rather strong dependence on 
the mode index n. Namely, the larger the number n the smaller the corresponding mode scattering 
lengths and as a consequence, the stronger is the scattering of this mode into the others. This strong 
dependence is due to squared transverse wave number = im/ d m the numerator and lengthwise 
wave number kn in the denominator of Eqs. (jH]) and (H5|) . Evidently, with an increase of the mode 
index n the value of fc„ decreases. An additional dependence appears because of the roughness power 
spectrum Vr(/c„=p/c„'). Since the binary correlator >V(x) of random surfaces is a decreasing function of 
its Fourier transform W{kn^kn') increases with n (note that it is constant for the delta-correlated 
roughness only). Therefore, all the factors contribute in the same direction for the dependence of Ln 
and on the mode index n. As a result, we arrive at the following hierarchy of mode scattering 
lengths: 



Li!f < LSi^i, < ... < < (46) 

The smallest mode attenuation lengths L^^j and belong to the highest channel with the mode 

index n = N^, while the largest scattering lengths and L^^ correspond to the lowest channel 
with n = 1. Note that similar hierarchy was also found in Refs. [Til Il2l S] in the model of quasi-one- 
dimensional systems with a stratified disorder. 

As is known, the quasi-one-dimensional systems with isotropic volume disorder reveal three 
typical transport regimes, the regimes of a ballistic, diffusive (metallic), and localized transport. In 
contrast to this conventional picture, in Refs. jlH] it was shown that in the case of surface disorder a 
very important phenomenon of the coexistence of ballistic, diffusive, and localized transport emerges. 
This happens due to hierarchy fH6|) of scattering lengths, and even in the absence of correlations in 
C,{x), Specifically, while the lowest modes can be in the ballistic regime, the intermediate and highest 
modes can exhibit the diffusive and localized behavior, respectively. This effect seems to be generic 
for the transport through the waveguides with random surfaces. 

One can see now that unlike the single-mode case, the concept of one-parameter scaling is no 
more valid for the transport in multi-mode surface-disordered systems. There are two points that 
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should be stressed in this respect. On the one hand, the average partial transmittances (T„) entering 
the equation ( H3|) . are very different for different conducting channels. On the other hand, and what is 
even more important, all propagating modes turn out to be mixed due to the inter- mode transitions. 
Therefore, the transmittance (T„) of a given n-th mode depends on the scattering into all modes, and 
the total average transmittance fj43l) is determined by the whole set of attenuation lengths f l4^ and 
(USD with l^n^Nd. 

Summarizing our brief discussion, it becomes clear that in quasi-one-dimensional guiding struc- 
tures with delta-correlated or Gaussian correlations in surface disorder, the crossover from the ballistic 
to localized transport is realized through the successive localization of highest propagating modes. 
Otherwise, if we start from the localized regime, the crossover to the ballistic transport is realized via 
the successive opening (delocalization) of lowest conducting channels. 

From this analysis one can conclude that for multi-mode structures with surface disorder the 
role of specific long-range correlations is much more sophisticated in comparison with single-mode 
waveguides. First, such correlations should result in the suppression of the interaction between differ- 
ent propagating modes. This non-trivial fact turns out to be crucial for the reduction of a system of 
mixed channels with quasi-one- dimensional transport to the subset of independent waveguide modes 
with a purely ID transport. Second, the same correlations can provide a complete transparency of 
each independent channel, similar to what happens in strictly ID geometry. 

To demonstrate these effects, let us take a random surface profile ^(x) of the form, 

i{x) = ^= / dx' Z{x - x') — , (47) 



with the slowly decaying (on average) binary correlator. It results in the "window function" for the 
roughness power spectrum. 



, , sinffcfX) 

Wix) = y \ (48a) 
W{k,) = ^e{k,-\k,\), k,>0. (48b) 

From equations (jH]) and (145|1 one can see that in the case of long-range correlations in a 
disordered surface (jlZj), the number of modes into which a given nth mode is scattered, i.e. the actual 
number of summands in the equations fl44p and 0451) . is entirely determined by the width k^ of the 
rectangular power spectrum fHHj) . It is clear that if the distance \kn — fcn±i| between neighboring 
quantum values of fc„ is larger than the correlation width kc, 



\kn kn±l \>kc, (49) 

then all inter-mode transitions (between different propagating modes) are forbidden. As a conse- 
quence, the sum over n' in the expression f H^ for the inverse forward scattering length contains only 
diagonal term with n' = n describing a direct intra-mode scattering inside the channels. Moreover, 
each term in the sum of Eq. f l45p for the inverse backscattering length is equal to zero. As a result, 
the following interesting phenomena arise. 

(i) All high propagating modes with indices n that satisfy the condition (H9l) . turn out to be 
independent of the others in spite of the interaction with rough surface. Therefore, they form a subset 
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( f) 

of ID non-interacting conducting channels with finite length of forward scattering L„ and infinite 

1 na^ {7in/d)^ 1 



backscattering length Ln^ 



0. (50) 



(ii) As is well known from the standard theory of one- dimensional localization (see, e.g., 
Refs. [201 Ell |22l |l2l HH |57]), the transport through any ID disordered system is determined only 
by the backscattering length Ln^ that in our consideration equals the localization length, and does 
not depend on the forward-scattering length L„ . Since the former diverges for every independent 



channel in line with the expression fl50|) . all of them are completely transparent because they exhibit 
the ballistic transport with the partial average transmittance (T„) = 1. This means that according 
to the Landauer's formula fH51) . the transmittance of the subset of such independent ballistic modes 
is simply equal to their total number. 

(iii) As for low propagating modes with the indices n contradicting to the condition (l49l) . they 
remain to be mixed by surface scattering because the roughness power spectrum fHSj) is non-zero for 
them, W{kn — kn') = n/kc- These mixed modes have finite forward and backscattering lengths and, 
consequently, stay in the diffusive or localized transport regime for large enough waveguide length L. 
As a result, they are non-transparent and do not contribute to the total transmittance (T). Therefore, 
the latter is equal to the number of independent ballistic modes. 

Note that the distance — kn±i \ between neighboring wave numbers kn and kn±i increases as 
the mode index n increases. Therefore, the inequality (H^ restricts the mode index n from below. 
That is why, in contrast with the conventional situation associated with the hierarchy of mode scat- 
tering lengths P6|) . the low modes are mixed and non-transparent, while high propagating modes are 
independent and ballistic. Because of a sharp behavior of the roughness power spectrum (jUj), the 
transition from mixed to independent modes is also sharp. 

More analytical results can be obtained for waveguides with large number of conducting channels 
Nfi, if the quantum numbers n of independent ballistic modes are also large: 

Nd = [kd/n] ^ kd/n ^ n > 1. (51) 

In this case the inequality (149|) is reduced to the requirement \dkn/dn\ > kc which can be rewritten 
in the following explicit form. 



n> K 



mix 



{kd/-K) 



(52) 



/I + (M/vr)- 

We recall that square brackets stand for the integer part of the inner expression. 

The condition fl52l) determines the total number Nmix of mixed non-transparent modes, the total 
number N^ai = — Nmix of independent ballistic modes, and the critical value of the mode index 
n that divides these two groups. All propagating modes with n > N^ix are independent and fully 
transparent, otherwise, they are mixed for n ^ N^ix and characterized by finite scattering lengths 
lI-P and Ln '' ■ Therefore, the total average transmittance fH51) of the multi-mode structure is. 



(T) = [kd/ir] - [kd/nad , = V 1 + (kcd/ir)-^ . (53) 

The numbers Nmix and N^ai of mixed non-transparent and independent ballistic modes are 
governed by two parameters, the mode parameter kd/n and the dimensionless correlation parameter 
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kcd/j:. In the case of "weak" correlations when kcd/j^ ^ 1, the number of mixed modes Nmix is of 
the order of N^, 



^ ' mix ~ 



kd\ 1 / kd\ ( krd 



TT 



TT / 



TT 



for kcd/rc ^ 1. 



(54) 



Consequently, in this case the number of balhstic modes N^ai is small, or there are no such modes 
at all. If the parameter kcd/n tends to infinity, kcd/n — t- oo, the rough surface profile becomes 
white-noise-like and, naturally, N^ix ~^ 

The most appropriate case is when a random surface profile is strongly correlated so that the 
correlation parameter is small, kcd/n <^ 1. Then the number of mixed non-transparent modes Nmix 
is much less than the total number of propagating modes Nd, 



kd 

TT 



kcd 



TT 



< Nd for kcd/TT < 1. 



(55) 



Therefore, the number of independent modes Nf,ai is large. When the correlation parameter k^d/TT 
decreases and becomes anomalously small, kcd/TT < (kd/TT)"^ ^ 1, the number of mixed modes N^ix 
vanishes and all modes become independent and perfectly transparent. Evidently, if the correlation 
parameter kcd/n vanishes, kcd/TT — )■ 0, the roughness power spectrum fHHl) . becomes delta-function-like 
and, as a consequence, Nmix = 0. In this case the correlated disorder results in a perfect transmission 
of waves. 

Finally, let us briefly discuss the expression (153|) for the total average transmittance. In Figure [3] 
an unusual non-monotonic stepwise dependence of (T) on the mode parameter kd/TT is shown, that 
is governed by the width kc of the rectangular power spectrum 




Figure 3: Step-wise transmittance fl53|) of surface-disordered guiding system versus the mode pa- 
rameter kd/TT. In the figure, the value of the normalized correlation parameter kcd/TT = 0.32 (after 

ID)- 



Let us discuss this Figure. Within the region where kd/TT < ac, the inter-mode transitions 
caused by specific surface correlations, are forbidden for all conducting channels. Therefore, all 
propagating modes are independent and ballistic, and the second term (the number of mixed modes) 
in the expression (153!) for the total transmittance is equal to zero. Here the transmittance exhibits 
a ballistic stepwise increase with an increase of the parameter kd/iT. Each step "up" arises for an 
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integer value of the mode parameter kd/n, when a new conducting channel emerges. Such a stepwise 
increase of the total transmittance is similar to that known to occur in quasi-one-dimensional ballistic 
non- disordered structures (see, e.g., |66]). 

Otherwise, when kd/n ^ Oc, in addition to the standard steps "up" originated from the first 
term in equation f lS5]) . there are also the steps "down" associated with the second term. These 
steps "down" are provided by the correlated surface scattering and arise when a successive low mode 
abruptly becomes mixed and non-transparent. In other words, the positions of the nth step "down" 
are at the transparency edge point kd/n = nac where the nth conducting channels closes. Specifically, 
the first step "down" occurs at the total transparency edge kdjn = ac, where the first mode is closed. 
This transparency edge separates the region of complete transparency from that where lower modes 
are mixed and non-transparent. The second step "down" is due to the particular transparency edge 
kd/n = 2q!c of the second mode, etc. Since the values of the ratio kd/nac are determined by the 
correlation parameter kc, the positions of steps "down", in general, do not coincide with those of steps 
"up" . The situation may also occur when the steps "up" and "down" cancel each other within some 
interval of the mode parameter kd/n. The interplay between steps "up" and "down" results in a new 
kind of stepwise non-monotonic dependence of the total quasi-one- dimensional transmittance. The 
experimental observation of the discussed non-conventional dependence seems to be highly interesting. 

F.M.I acknowledges the support by CONACyT grant No. 80715. 
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